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Abstract 

We explicitly construct a large class of unitary transformations that allow to perform 
the ideal estimation of the phase-shift on a single-mode radiation field. The ideal 
^ ^ phase distribution is obtained by heterodyne detection on two radiation modes after 

■ the interaction. 

o 

m 

^ ■ The quantum estimation of an unknown phase shift — the so called quantum 

p ■ phase measurement — is the essential problem of high sensitive interferometry, 

and has received much attention in quantum optics [1]. For a single-mode 
C \ electromagnetic field, the measurement cannot be achieved exactly, even in 

^ \ principle, due to the lack of a unique self-adjoint operator [2]. In fact, the ab- 

CT' sence of a proper self-adjoint operator is mainly due to the semi-boundedness 

of the spectrum of the number operator [3,4], which is canonically conjugated 
^ , to the phase in the sense of a Fourier-transform pair [5]. 



This observation opened the route for an exact phase measurement in terms of 
two-mode fields, where a phase-difference operator becomes conjugated to an 
unbounded number- difference operator [6]. In fact, a concrete experimental 
setup using unconventional heterodyne detection has been suggested [7] for 
this kind of measurement. However in the single- mode case, no feasible scheme 
that can provide the optimal phase measurement has been devised yet. 

The most general and concrete approach to the problem of the phase measure- 
ment is quantum estimation theory [8] , a framework that has become popular 
only in the last ten years in the field of quantum information. Quantum es- 
timation theory provides a more general description of quantum statistics in 
terms of POVM's (positive operator-valued measures) and gives the theoreti- 
cal definition of an optimized phase measurement. The most powerful method 
for deriving the optimal phase measurement was given by Holevo [9] in the 
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covariant case. In this way the optimal POVM for phase estimation has been 
derived for a single-mode field. More generally, the problem of estimating the 
phase shift has been addressed in Ref. [10] for any degenerate shift operator 
with discrete spectrum, either bounded, bounded from below, or unbounded, 
extending the Holevo method for the covariant estimation problem. 

As already stated, quantum estimation theory provides the optimal POVM 
for the phase measurement. This writes in terms of projectors on Susskind- 
Glogower states [11] 



c^M<^) = ^|e^")(e'1, (1) 

where |e*'^) = EJ^o Notice that the states |e*'^) are not normahzable, 
neither orthogonal, however they provide a resolution of the identity, and thus 
guarantee the completeness of the POVM, namely 



rd^{^)^I (2) 
io 

For a system in state p, the POVM in Eq. (1) gives the ideal phase distribution 
p(</7) according to Born's rule 



p{^)=T,[d^{^)p\ = ^{e^W^). (3) 

In this Letter we will explicitly construct some unitary transformations that 
allows to perform the ideal phase measurement described by the POVM in 
Eq. (1). First, we will introduce an isometry V which enlarges the Hilbert 
space of the system (say Tia for mode a) to the tensor product TY^ ® 'Hh for 
two modes a and h. Then, we will prove that the exact measurement of the 
complex photocurrent Z = a — l)^ provides through its marginal distribution 
the ideal probability density p(<^) of Eq. (3). Finally, we will construct a large 
class of unitary operators on Ha ® Hb <Si He, where He denotes the Hilbert 
space of an ancillary arbitrary system, such that the isometry V is realized 
with unit probability. 

We start by introducing the eigenstates of the heterodyne photocurrent Z — 
a-b'' [12,13,7] 



Z\D{z)))al> = z\D{z)))al> , (4) 

where D{z) = exp(2;a''^ — z*a) denotes the displacement operator. Here and in 
the following we use the notation [14] for bipartite pure states on Ha ® Hb 
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\^))ab= E Anm\n)a®\m), = A® h\I))^, = h® A^\I))^, , (5) 



n,m=0 



where denotes the transposed operator with respect to some pre-chosen 
orthonormal basis. The states \D{zyf)ab are orthogonal in Dirac sense over the 
complex plane, namely 



ab{{D{a)\D{(3)))ab = T^5^'^\a -(3) = 7r5(Re a -Re (3) 5(Im a - Im /3) .(6) 
They also provides a basis for Ha ® Ti-b as follows 



— \D{z)))abab{{D{z)\^Ia®Ib. (7) 



C TT 




The measurement of the complex photocurrent Z can be performed through 
unconventional heterodyne detection [13] with both the signal a and the image- 
band b non-vacuum (in usual heterodyne detection the image-band mode is 

in the vacuum, thus providing the well-known coherent-state POVM). The 
measurement of Z is also equivalent to two separate homodync measurements 
on modes -^(6 ± a). In fact, consider the 50/50 beam splitter operator R = 
exp[^(a^6 — a6^)] that realizes the unitary transformation 



(8) 



Upon denoting with | x)a and | y)b the eigenstates of the quadratures Xa — 
{a + a))/\/2 and = (ift'*' — ib)/\/2, one has the following identity [15] 

R{\x)aa{x\®\y)bb{y\)R^ = \D{x + zy)))abab{{Dix + zy)\ . (9) 

Notice also that this kind of measurement is performed in the teleportation 
protocol for continuous variable of Braunstein-Kimble scheme [16,17]. 

We can now write the isometry V such that the transformation 

Tip) = VpV^ (10) 

maps the state of the system p e Tia to a two-mode state in Tia ® Ti-b- The 
operator V has the form 

y-^l d^af{\a\)\D{a)))aba{e'^'^''\ ■ (11) 
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By choosing f{t) as an arbitrary function satisfying the condition 



f + OO 1 

/ dtt\f{t)\' = -, (12) 

Jo TT 

it follows that V is an isometry, namely V^V — la- 

It is easy to check that the transformation (10) has the following covariance 
symmetry 

r(e^^"^Ve-*^"^") = e^^''^'' ^ e-''^^^^T{p) e"'^"^" ^ e'^''''^ . (13) 

We can now evaluate the probabihty density of getting outcome 2; e C through 
the measurement of the photocurrent Z. One has 



TT 

-\\!m\\k^'''W)a. V^^rgz (14) 

From condition in Eq. (12), it follows that the marginal distribution on the 
statistical variable (p — arg z corresponds to the ideal distribution of the phase 
for any input state p. 

In the following we explicitly construct some unitary realizations of the map 
in Eq. (10). We start by defining the operators in C{7ia ® Hb) 



V=V{Ia®b{x\), -{Ia^\x)b)V^ , (15) 

where \x) is an arbitrary normalized state in Tib, and the tensor notation 
^a'^bixl represents a linear operator from Tia®Ti.b to Tia (the bra ^(xl can be 
regarded as a linear functional from Tib to C). Similarly, Ia®\x)h represents an 
operator from Tia to Tia ® V-b- Notice that both VV"^ and VW are projectors, 
namely 



vvWv^ = vv^ , vVvV = . (16) 

Upon introducing an arbitrary Hilbert space Tic (also finite dimensional), we 
construct the following operator 



U^V® WW^ -V^ ® W^W + V^V) ®W^ + {I- VV^) (g) W ,(17) 
where is a linear operator in Tic- Under the conditions 



4 



(18) 



one can easily check that WW^ and W^W are projectors orthogonal each 
other, and U is unitary. Consider now the transformation of the system pre- 
pared in a state p®(7®^ which has been evolved through the unitary U and 
traced over the ancillary space TYc- One has 

Trc[C/(p (g) cr ® n)U^] = V{p ® (t)V'^T^[WW^ n] 

+V{p 0a)(I- VV^)Tr[W^iJ,] + 0(p ® a)VTr[WWii] 

-y^(p(8) (7)(7 - V^V)Tr[Wi4 -(I- V^V)(p® a)VTr[W^ ii] 

+(I - yV)(p 0a){I- V''V)Ti[WW'' ii] + VV'')(p a)V''Tr[Wij] 

+(I -VV^)(p®a)(I -VV^)Tr[WWli] . (19) 

The map in Eq. (10) can then be achieved by the unitary transformation U, 
by taking 

^=\x)bb{x\, (20) 
and p, such that 

Tr[iyiyV] = 1 , Tr[Wp] = Tr[W^p] = . (21) 

We summarize the conditions on the measurement scheme: f{t), \x),W,p in 
Eqs. (11), (15), (17), (19) are arbitrary, provided that conditions (12), (18), 
(21) are satisfied. 

As an example, consider the case in which the space He pertains to a radiation 
mode c. One can take 



W=J2\2n)ee{2n + l\, p=\0),,{0\, (22) 

n=0 

thus obtaining 



T{p) = Tre[U{p®\x)bb{x\ ® |0)cc(0|)t/1 . (23) 

We notice that the "pseudo-spin" operator W in Eq. (22) has been introduced 
also in Refs. [18,19] in the context of Bell's inequalities for continuous variable. 

The result in Eq. (23) is similarly obtained for a qubit system in Tic, with 

|0)cc(l|. 
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In conclusion, a large class of unitary realizations of the ideal phase mea- 
surement of a single-mode radiation field has been presented. These unitary 
evolutions act on the Hilbert space Tia ^V-b® Ti-c, where Ha and Hb are re- 
ferred to radiation modes, and He pertains to an arbitrary system. By suitably 
preparing the state of the systems in Hb and He, the ideal phase distribution 
for the input state p G Ha is obtained through heterodyne detection performed 
after the interaction on modes a and b. 
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